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Abstract 

We investigate finite energy solutions of Yang-Mills-Chern-Simons systems in odd spacetime dimen- 
sions, d = 2n + 1, with n > 2. This can be carried out systematically for all n, but the cases n = 3, 4 
corresponding to a 7, 8 dimensional spacetime are treated concretely. These are static and spherically 
symmetric configurations, defined in a flat Minkowski background. The value of the electric charge is 
fixed by the Chern- Simons coupling constant. 



1 Introduction 



Yang-Mills (YM) fields in d = D+l dimensional spacetime have interesting properties, especially in the static 
limit where the solutions to the equations of motion describe finite energy topologically stable solutions. 
Famously, in the case D = 3 these are the 't Hooft-Polyakov monopoles [TJ [3] on R 3 whose topological 



(f) \ charge is the magnetic flux. In fact, monopoles exist also on all R D |3]. Monopoles are supported by Yang- 



Mills-Higgs (YMH) systems, where the Higgs field (which is always an iso-Z?— vector) defines the topology 
completely by virtue of the requirement of finiteness of energy [H [5] . 

However, in d = 3 + 1 dimensional spacetime, there exists also the Julia-Zee dyon [6] which in addition to 
magnetic flux also possesses an electric flux. Unlike the magnetic charge, the electric charge of a dyon, while 
it is a global charge, is not a topological charge. The topological charge appearing in the Julia-Zee (JZ) 
dyon is the magnetic monopole charge given by the spacelike component Ai of the gauge connection (and the 
Higgs field $), while the electric flux is given by the timelike component Aq (and <£>). By construction, Aq 
and $ are proportional in the absence of a Higgs potential. Now the existence of the Julia-Zee (JZ) dyon is 
predicated on the presence of a monopole, as well as the availability of the 't Hooft electromagnetic tensor. 
This mechanism is restricted to 3 + 1 dimensional spacetime. Although monopoles can be constructed in all 
dimensions, the definition of a 't Hooft tensor in higher dimensional space is problematic, as discussed in 
[3]- So it is not open to us to exploit the I?— dimensional monopoles for the purpose of constructing dyons 
in higher dimensions. 

To the best of our knowledge, the only dyon in higher dimensions is the dyonic instanton of Lambert and 
Tong [14 j (LT) in 4 + 1 dimensional spacetime. The topological charge of the LT dyon is the usual instanton 
number on the 1R subspace, described by the spacelike component Aj of the gauge connection. The electic 
flux is again given by the timelike component Aq and the Higgs field <&, as for the JZ dyon. Unlike in the 
case of the latter however, Aq of the LT dyon is actually equal to $ and not just proportional to it. The 
dyonic instanton describes a global electric flux in addition to the Pontryagin charge (of the second Chern 
class), the latter being the topological charge analogous to the magnetic charge of the Julia-Zee dyon. 

Our aim in the present work is to construct a new type of instantonic dyons in (higher) odd dimensional 
spacetimes. For these solutions, the spacelike sector is stabilised by a topological charge, and which support 
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also an electric flux. As in the case of the LT dyonic instanton, where the "magnetic charge" is the 2— nd 
Chern-Pontryagin charge, the topological charge here is the n— th Chern-Pontryagin charge on ~R D , D = In. 
The latter is given by the spacelike, or "magnetic" , component Ai of the gauge connection. Our solutions 
share this property with the LT dyonic instanton. Naturally, in our case too, the electric flux associated 
with the solutions proposed is given by the timelike, or "electric", component A . 

The models proposed here differ from those of JZ [B] and of LT [T3] in two respects. In common with 
these, they feature non-Abelian Yang-Mills fields, but in contrast they contain no Higgs field. The Higgs 
field is employed to support the magnetic (topological) charge of the JZ dyon, and in the case of the LT 
dyonic instanton, it supports the component Aq . The role corresponding to the Higgs (kinetic) term is here 
played by the non-Abelian Chern-Simons (CS) density. Employing a CS term in the action for the purpose 
of supporting Aq is standard, both for Abelian [TSJ [TB] and non-Abelian [T7] systems in 2 + 1 dimensions, as 
well as in higher dimensions j!8j . 

It is important to emphasise the difference in using the instanton number as the topological charge, 
rather than the monopole charge. The connection A, of the instanton behaves asymptotically as a pure- 
gauge, in contrast to that of monopole which behaves as half pure-gauge. In other words the latter is a 
Dirac-Yang [71 [HI IS] monopole decaying as r _1 , more slowly than the instanton. As a consequence, the 
energy density functionals of monopoles in higher (than three spacelike) dimensions cannot involve the usual 
quadratic Yang-Mills term whose energy integral diverges. This is not the case for the faster (pure-gauge) 
decay of the instanton, which enables the retention of a quadratic YM density in all dimensions, with 
converging "energy integral". From a physical standpoint, this alone could be considered a motivation for 
employing instanton numbers in preference to monopole charges as topologiccal charges in higher dimensions. 

The models we will introduce will feature the Yang-Mills sector, consisting of some or all possible terms of 
the Yang-Mills hierarchy jTTjrTS] (to be introduced below), and the Chern-Simons density in the appropriate 
dimension. The topological charge is the instanton number of the given YM system on the "spacelike" 
subspace H D (or R 2 "), which always includes the usual (quadratic) Yang- Mills density, as well as at least 
one other member of the YM hierarchy which is of sufficiently high order in the YM curvature to enable the 
Derrick scaling requirement for the convergence of the "energy integral" to be satisfied. The application of 
the Derrick scaling requirement Yang-Mills sector (on a Euclidean space) is rigorous. 

It is important to distinguish the status of the Derrick scaling requirement in the case of dyons where 
the electric YM connection Aq is introduced in the covariant derivative of the Higgs field as is the case for 
the JZ and LT dyons, and in our case when it is introduced via a CS term in the action as was done in the 
case of the gravitating YMCS system in 4 + 1 and higher dimensions [18j Qjj] . (Note also the Abelian [TH [16] 
and non-Abelian [17] CS-Higgs vortices in 2 + 1 dimensions, where Aq enters both in the Higgs covariant 
derivative and the the CS density.) The status of Derrick scaling is subtle when there is a CS term in the 
action, and the non-Abelian case is less transparent than the Abelian. In the Abelian case, the electric gauge 
connection Aq can be found explicitly by solving the Gauss Law equation, and substituting it into the static 
energy density functional exposes the (Derrick) scaling properties. In the non-Abelian cases [ITJ [T8] [19] by 
contrast, it is not practicable to solve the Gauss Law equation to yield the (non-Abelian) electric gauge 
connection Aq explicitly. In those cases the solution of the Gauss Law equation is implicit in the numerical 
process, resulting nonetheless in the required scaling contribution of the CS term appearing in the action, 
now in the static energy density. Thus the CS density acts as a higher order term enabling the Derrick 
scaling requirement for the finiteness of energy to be satisfied. 

Our formulation is in principle for all spacetime dimensions d = 2n + 1, n > 3. We have considered 
spherically symmetric static solutions and the asymptotic analysis for these is given for arbitrary D = In. 
Explicit, numerical, constructions however are restricted to the n — 3, 4 cases. 

In Section 2 the models are introduced, and the imposition of spherical symmetry as well as the resulting 
field equations are given in Section 3. The solutions are presented in Section 4, which includes the asymp- 
totic analysis in the general case, as well as the concrete numerical construction for the dyons in 7 and 9 
dimensional spacetimes. In addition, we supply three Appendices. Appendix A presents a simplification 
of the (dynamical) Chern-Simons density in the static limit. This simplifies the algebraic calculations con- 
siderably, in particular enabling the formulation of the arbitrary dimensional case. Appendix B gives the 
spherically symmetric Ansatz for the full SO(D + 2) Yang-Mills system. Appendix C gives an attempt at 
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generalising the Lambert- Tong 14! dyonic instanton in.<f = 4 + ltod = 4p+l. While this attempt is not 
completely satisfactory, the review of that material exposes the essential difference of our instantonic dyons 
in 7, 9, . . . dimensions and the LT dyonic instanton in 5 dimensions. 



2 The models 

In selecting the models in d = D + 1, with D = 2n, we have invoked two criteria, that of topological 
stability of the "magnetic" sector, and, that of the existence of a global electric charge. Our first criterion 
is that the spacelike components of the non-Abelian gauge fields describe a topologically stable (static) field 
configuration charcterised by an instanton charge. (The usual action integral here is interpreted as the 
"energy integral".) For these configurations to have finite energy, it is necessary to employ at least two 
appropriately scaling members of the Yang- Mills hierarchy [TTJ [T2J [3] on H D = R 2 " 



in terms of the "magnetic" field Aj. Later, when we introduce the "electric" field Aq, this definition will be 
retained, but then on d = D + 1 dimensional spacetime with Minkowskian signature. In (|2.1j) . F{2p) is the 
totally antisymmetrised p— fold product of the YM curvature 2— form F(2). We shall refer to the system 
(|2.1j) as a p-YM system, the 1-YM system being the usual YM density. Note that the p-YM density scales 
as L~ ip . In any given spacetime dimension d, the constraint of antisymmetry of F(2p) requires that the 
highest order curvature term F(2P) is that with d = 2P, which scales as L~ iP . 

Subject to this constraint, the most general YM system is the superposition of terms (|2.1|) 

£ Y M = £^Tr{w}, (2.2) 

the p = 1 term being the usual YM system. Finiteness of the "energy" requires that Derrick scaling be 
satisfied, and given that the lowest order curvature term with p = 1 scales as L~ 4 and the highest with 
p = P as L~ 2 ^ D+1 \ this is (more than) sufficient for Derrick scaling to be satisfied on H D . A subsystem of 
the superposition (|2.2[) will be adopted as the Yang-Mills sector of our models. 

Our choice of the Yang-Mills sector is made such that the solutions on H D , i.e. the "magnetic" field 
configurations, have finite "energy" and be topologically stable, as is the case for the Julia-Zee dyon and the 
dyonic instanton. The topological stability stems from the following sets of inequalities 

Tr(F(2 Pl )~(*F(2p 2 ))(2p 1 )\ >0, (2.3) 



where *F(2p 2 ) is the Hodge dual of F(2p 2 ), which is of course a 2pi-form, provided that the (pi,p 2 ) pah 
is a partition of 2{p\ + p 2 ) = D = 2n. This means that topological stability constrains the highest order 

curvature term in (|2.2[) . Tr^F(2P) 2 j> be the P = D term, rather than P = D + 1 in Minkowski space. 

It follows that for any D = 2n, one has the "energy" lower bound 

/•(Pi<P2) def. r ( Pl ) r ( P2 ) ^ r (O A\ 

for any partition^ n = p\ +p 2 ■ Here, Cua is the n— th Chern-Pontryagin density, and (ri , t 2 ) are dimensionful 
(coupling) constants. 

The criterion of topological stability (12.41) requires the presence of at least two YM terms. Thus, for 
simplicity we shall restrict our definitions to systems consisting of exactly two p-YM terms. Moreover for 



1 When pi = p2 = p, with D = 4p, the topological inequalities (I -'.'111 and (l-'.4t) are saturated and there result self-dual 
BPST | 10 | configurations |lll . the p = 1 case being the usual BPST instanton. We have eschewed this choice here since no 
solutions exist for such systems when the Chern-Simons terms are introduced. A discussion of this will be given in Section 5. 
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"physical" reasons, we will fix p\ to p\ = 1, so as to retain the usual (quadratic) 1-YM term. This in turn 
fixes fix f>2 to P2 = n— 1, finally fixing the YM system to the superposition of the 1-YM and the (n — 1)-YM 
terms. 

Finally, the existence of 'instantons' of these systems requires that the YM connection takes its values in the 
chiral representation of SO(D) = SO(2n), such that the gauge group G must contain the subgroup SO(D). 
This completes the definition of the Yang-Mills sector. 

Concerning the introduction of the Chern-Simons (CS) term (which fulfills our second criterion, namely 
that of supporting an electric field,) this is uniquely fixed by the dimension D + 1 of the spacetime, and 
is accompanied with the introduction of the "electric" connection Aq to the whole system. The smallest 
simple gauge group that supports a nonvanishing CS density is G = SO(D + 2), which finally fixes our 
choice of gauge group. In practice the gauge group will be truncated to G = SO(D) x SO{2) because in 
the spherically symmetric case studied here, it transpires from the asyptotic analysis in Section 4 that only 
G = SO{D) x SO{2) solutions can be found. (We expect that in the presence of a negative cosmological 
constant there exist full SO(D + 2) solutions, as was found in [T8 for d = 5. It is for this reason that the 
full SO(D + 2) spherically symmetric Ansatz is given in Appendix B.) 

As in all dyonic systems, the topological stability of the purely "magnetic" sector does not any more 
guarantee the existence of finite energy solutions. (Indeed analytic proofs of existence for dyons are prob- 
lematic and their existence is usually established numerically, as for the JZ dyon (in the presence of a Higgs 
potential), or by explicit construction in closed form as in the case of the LT dyonic instanton.) 

Therefore the Lagrangians we adopt are 

£yMCS = £YM + ) ( 2 - 6 ) 

in which £ym is given by (|2.5|) and C^g is the CS density in d = 2n + 1 dimensional spacetime. 

The definition of the Chern-Simons densities is standard, if complicated, but here we are interested only 
in static field configurations, in which case it simplifies considerably. The relevant formulas are presented in 
the Appendix, and here we simply state the simplified definition of the static CS densities. With dtA^ = 0, 
one can show that, up to a total divergence term (which we ignore here since we are only interested in the 
Euler-Lagrange equations) , the effective arbitarary n CS Lagrangian reduces to the effective density 

4" s } = (n + 1) ^w^.^-^Tr j Aq Fisu F.^^ | . (2 . 7) 

To our knowledge, this is a new result, and its derivation is given in Appendix A. 

Finally, we define the scalar valued global charges of our solutions. The global charges of a dyon are the 
magnetic and the electric fluxes. The "magnetic" flux on R 2n here is the n-th Chern-Pontryagin charge C^ n \ 
appearing e.g. in (|2.4p and (|2.3[) . The definition is familiar, in terms of the static "magnetic" curvature Fij. 

Next, definition of the electric flux. As always, the appearance of a Chern-Simons term in the Lagrangian 
gives rise to a nontrivial timelike, or "electric" , component of the gauge field. The electric flux here 

dS> Ei , (2.8) 

(Ei = Fio) is a non-Abelian quantity. However, this flux takes its values always in one single element of the 
SO(D + 2) algebra, and can therefore be interpreted as the global electric charge. The reason is that at 
infinity the field is radially symmetric, and subject to spherical symmetry the the solutions are restricted to 
the SO(D) x SO (2) subalgebra of SO(D + 2), and, A takes its values along the 50(2) subalgebra. This 
fact will become clear from the asymptotic analysis below (see (|4.26jl ). 

An alternative definition for the electric flux can be given such that it is expressed as a scalar valued 
global charge, ab initio. This is in the same spirit as in [T3], the only difference here being that we do not 
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have a Higgs field, instead of which we employ the non-Abelian gauge connection Aq 



<l 



dcf. 



J dSiTrlAoEiX = J dSi Tr A D L A Q = i ^ dS^Trj^j^ ^ d n (23) ±. A 2 Q , (2.9) 



where dflm) is the angular volume element in R . Both definitions ()2.8[) and (|2.9|) give, up to normalisation, 
the same result. 

Before proceeding to the consideration of specific models considered, we return to our statement in 
footnote 1 , namely that the choice of YM sector £ym in (|2.6I) , given by p\ = P2 = in (|2.4[) , is excluded 
in the present work. There exist no finite energy solutions to those systems, but as seen in the D = A 
(n = 2) case in 5-dimensional spacetime such solutions can be constructed when a suitable scalar field is 
introduced [20] . Whether such solutions can be constructed in higher dimesnional spacetimes with n > 3 by 
the introduction of a scalar field is an open question. 

In our concrete numerical constructions, we will consider the following two actions, in d = 6 + 1 and 
d = 8 + 1 spacetime dimensions, n = 3 and n = 4, with D = 2n. 



S 7 = 
S 9 = 



/ d 7 Xy/—g 

J M 



d xy/^g 



M 



2 • 2! 

n 
2 • 2! 



2-4! 
2 • 6! 



9^ r( 4 ) 



M 



cs 



(2.10) 



(2.11) 



where = cJ^A^ — + [A M ,.A„] is the gauge field strength tensor, and k is the CS coupling constant. 



3 The Ansatz and field equations 

In Appendix B, the static SO(D + 2) YM system on M D (D = 2n) is subjected to spherical symmetry. In 
the present work, we use a particular truncation of this Ansatz, namely one where the functions and x 
(which define the magnetic and electric potentials, respectively) are replaced by 

<£= (to, 0,0) and x = (V,0,0). 

This is because subject to spherical symmetry, we could construct only SO{D) x 50(2) solutions, and not 
the ones for the full SO(D + 2) gauge group. We have presented the general case in Appendix B since we 
expect that introduction of a negative cosmological constant (i.e. an Anti-de Sitter spacetime background) 
would enable the systematic extension of the solution presented here to the full SO(D + 2) case □• The 
SO(D) x SO (2) truncated version of the Ansatz employed in this work is 

tot V ) 1 X* 

A=-^ E«— dx 1 +V(r)T, D+1D+2 dt, with i,j = 1, . . . , D , (3.12) 

r r 

being the chiral representation matrices of SO(D), and Eu + i ! ij + 2 of the SO(2), subalgebras in SO(D+2), 
defined in Appendix B. r and t are the radial and time coordinate in the (D + l)-dimensional Minkowski 
space, while x l are the usual Cartesian coordinates, with x l x % = r 2 . it is convenient to use the following 
chiral Sigma matrices in (D + 2) dimensions, 

Yii = —Tii — iTi , Kfl+i = — S_d+i = iTn+x , = = I (3.13) 

in terms of the gamma matrices Tm — (r^Tm) in (D + 2) dimensions. 

2 In |18| (for D = 4), the presence of a negative cosmological constant resulted in the full SO(D + 2) solutions being realised 
for the gravitating YMCS system. 
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Subject to the spherical symmetry Ansatz (|3.12[) the p-YM ansatz density on H D reduces to the following 
one dimensional functional of r 



Tr{ F(2pf 



(D-l)l rD ^(w*-i 



(D-2p+l)\ 
+(2p-l)(£>-2p+l) 



2(p-l) 



(2p) V' 2 + 



(2p) 



(D - 2p) 



w 2 - 1 



(3.14) 



where the power of r D_1 in the volume element is included. 

The corresponding one dimensional residual density pertaining to (|2.7|) for arbitrary D — 2n can readily 
be calculated 



L 



(n) 

(CS) 



n(n + l)V(w 2 -l) 71 -^' 



(3.15) 



Substituting this Ansatz in ()2.6|) . we find the one-dimensional reduced Lagrangian 



'e// 



dtdr 



„D-1 



■iv- + (D-D 



2r 2 



(Z) — 2)(1 - w 2 
4^4 



(3 



1-w 2 



2(2p-l)! 
kD(D + 2)V(w 2 - l)^ D -^w' 



-V u + (2p-l)(D-2p+l) ( — 



W" (D — 2p)(l — w ) 



2pr 4 



where a prime denotes a derivative with respect to r. Here we have scaled out the factor t\ and denoted 
t = T2 uZzp)] i a ^ so t° simplify the relations we have denoted P2 = P henceforth. 
The corresponding energy of the solutions is given by the integral 



E = Vd-i 



dr 



„D-1 



„D-1 



V' 2 + (D-1) 



2r 2 



(D-2)(l -w 
4^" 



2\2 



(3.17) 



1 - w 1 



2\ 2(p-l) 



2(2p-l)! 

The resulting (ordinary differential) equations are 
(D-2p+l) 1 1-w 2 



V rW + (2p-l)(X>-2p+l) ( — + 



(£>-2p)(l- w 2 ) 2 



2pr 4 



w" 1 + T- 



l-^ 2 a(p _ x A (13 - 3)w> _ (D-2)w(w 2 -l) 
(D-l) (2p-2)r r 2 ' J r r 2 

(2p-l)! £>-l 1 r 2 ' \ + 



(3.18) 



D(Z) + 2)k dF(iti)^, ^ 1 2(p-l),w 2 -l 



(£> - 1) rfw 



(2p-l)(D-2p+l) 
2(p - l)r 4 



(rw'{{D - 4p + l)(w 2 - 1) + 2(p - l)rW) - (D - 2p)w(l - w 2 ) 2 ) = 0, 



and 



1/' 1 



.(1z^!)2(p-i) 



+ 2)F(u>) 



0. 



(3.19) 



(2p-l)P r 2 

the last one being the Gauss Law equation. The function F(w) in fl3.18[l - (|3.19[l has the following general 
expression in terms of the hypergeometric function 2 -F\: 



, 1 2- D 3 
= -(-1) =■ 2-Fi ( 2'^ — '2 ;ur 1 "' 



(3.20) 



As seen from the relation (|3 . 19[) . an important generic feature of YMCS models within the SO(D) x 50(2) 
truncation is the existence of a first integral for the electric potential V(r), 



V 



1 



P-kD{D + 2) 



r D-l -i , r ( w 2 -l \2(p-l) 



jri»=0, 



(3.21) 



G 



with P an integration constant, which, as we shall see, for regular solutions is fixed by k. 
One can note that equations (|3.18[) and (|3.21l) are invariant under the scaling 

r^Xr, w^w, V^V/X, t ->• A 4 ^ -1 V, and k^A^k, (3.22) 

with A an arbitrary positive parameter. Then, without any loss of generality, one can set in this way r or k 
to take an arbitrary nonzero value. 

Another symmetry of the equations of the model p,18|) . (|3.19p consists in simultaneously changing the 
sign of the CS coupling constant together with the electric or magnetic potential 

k — > —K, V — > —V, or k — > —K, w — > —w. (3.23) 

In what follows, we shall use this symmetry to study solutions with a positive k only. 



4 The solutions 

We start by presenting the expansion of the solutions at near the origin r = 0. The regularity of the gauge 
field implies w — > ±1 there. Then it follows that the parameter P in the first integral (|3.2ip is fixed by the 
value of the CS coupling constant. Technically, this results from the fact that the term P — nD(D + 2)F(±1) 
should vanish as r — > 0. Restricting without any loss of generality to w(0) = 1, one finds 

P = kD(D + 2)F(1). (4.24) 

Then one finds the following expansionsnear the origin: 

"M-i-*+°e>. ^>=' (£ff^%;V +0 ' f ' ) ' (4 - 25) 

The only free parameters here is b = —hw"(0). The coefficients of all higher order terms in the r — > 
expansion are fixed by b. 

The solutions have the following expansion^] as r — > oo: 

w(r) = -l--^ + 0(l/r D ), V(r)=V -^ + O(l/r D ). (4.26) 

In the above relations, J, Vo are parameters given by numerics which fix all higher order terms, while Q is 
a constant fixing the electric charge of the solutions, 

Q = Ki . 1} ^^>^. y (4 , 7) 

The solutions interpolating between the asymptotics (|4.25l) , (|4.26l) were constructed by using a standard 
Runge-Kutta ordinary differential equation solver. In this approach we evaluate the initial conditions at 
r = 10 -5 , for global tolerance 10 -12 , adjusting for shooting parameters and integrating towards r — > oo. 
The input parameters are r and k (note that only the ratio of these two constant is relevant). 

Our numerical analysis gives evidence for the existence of nontrivial solutions of the eqs. (13.181) . (|3.19[) for 
D = 6, 8, these cases being studied in a systematic way. However, solutions with similar properties should 
exist for all higher values of D = 2n. The profile of a typical solution for the case D — 6, p = 2 is given in 
Figure 1. For all solutions we have constructed, the gauge function w(r) has a monotonic behaviour, with 
a single node for an intermediate value of r. A monotonic behaviour has been also noticed for the electric 
potential V(r). For moderate values of the ratio T2/K, the energy density is mainly concentrated in a shall, 
its maximum corresponding to the position of the of the node of the magnetic gauge function. 



J Note that although ui(oo) = 1 is also allowed, we could not find such multinode solutions. 
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Figure 2. The total energy E, the shooting parameter b — — |iu"(0), and the asymptotic value Vo of the elec- 
tric potential are shown as a function of the parameter k (left) and as a function of the parameter T2 (right) for 
D = 6, 8 solutions, respectively. 

The dependence of solutions on the parameters k, r 2 is shown in Figure 2. Starting with the dependence 
on the parameter k, we notice that nontrivial solutions with finite mass are also found in the absence of a 
CS term, n = 0. However, these configurations have a vanishing electric potential and correspond to the 
instantons in [T3]. The total mass of solutions increases almost linearly with k, while the shooting parameter 
b = —w"(0)/2 decreases. Also, as one can see in Figure 2 (right), one can find finite energy solutions of the 
YMCS equations even in the absence of the £^ 1 ' term in (|2.5[) . i.e. for r 2 = 0. 



5 Summary and discussion 

We have formulated the construction of dyon solutions to a family of Yang-Mills-Chern-Simons (YMCS) 
systems in 2n + 1 (n > 3) dimensional flat spacetime. These are static spherically symmetric field config- 
urations carrying both magnetic and electric flux. The components of the non-Abelian gauge connection 
= (Ai,A ) are denned to be in the algebra of SO(D + 2), but subject to spherically symmetry the 
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solutions turn out to take their values in the subalgebra SO(D) x 50(2). 

Just as for the other known dyonic solutions in four [B] and five |14j dimensional spacetimes, the "mag- 
netic" flux is a topological charge. In four dimensions the topological charge of the Julia-Zee (JZ) dyon is 
the monopole charge on R 3 , which is a descendant 3 of the second Chern-Pontryagin (CP) charge. In five 
dimensions the topological charge of the Lambert- Tong (LT) dyonic instanton, it is the second CP charge 
itself on R 4 . In our case here, the topological charge is the n-th CP charge on R™, n > 3. In the above, we 
have referred to the topological charges as "magnetic" charges, in analogy with |Bj, and also because these 
are defined exclusively in terms of the static spacelike (or "magnetic") components of the connection Ai, 
with the timelike (or "electric" ) components Aq entering the definition of the "electric" flux. 

In this context the JZ dyon can be described as a monopolic dyon while the LT dyon and our solutions 
can be described as a instantonic dyons. The salient difference between the two types of topological charge 
is that the former decay like Dirac-Yang [9] monopoles and hence in all but 3 + 1 dimensions the energy 
integral of the usual (quadratic) Yang-Mills (YM) term diverges, while the latter decay (like instantons) 
faster and the "energy" integral of the usual YM term converges. This can be seen as a desirable 'physical' 
feature. 

Concerning the detailed choice of the models proposed, the gauge field system supporting the topological 
("magnetic") charge is directly specified like in the known dyonic models. In the case of the JZ dyon, 
this is the usual SU(2) Georgi-Glashow model 0. In the case of the LT dyon, this is the usual (quadratic) 
SO±(4) — SUl/r(2) Yang-Mills model. In our case, like LT, the sector supporting the topological charge is 
described by the non-Abelian gauge field only, and the natural system is then the superposition of members 
of the Yang-Mills hierarchy (|2.2I) . Now the number of terms in (|2.2I) grow with increasing D (or n) and it 
may be reasonable to retain the minimum number of these necessary for the Derrick scaling requirement to 
be satisfied. For this, one can select any pair of terms, both with even p, for which the topological inequality 
(|2.3|) can be satisfied. Of such pairs, it may be reasonable to privilege the pair pi = 1 , p2 = n — 1, if only 
for "physical" reasons favouring the presence of the usual (quadratic) YM term. This is not necessary, but 
it is what we have done and as it happens in the two examples in 7 and 9 dimensions studied quantitatively, 
these are the unique possibilities. 

Alternatively, one might have opted to retain only the p = ^ term in (12.21) . i.e. with t\ — in (I2.5[) . 
which on its own satisfies the Derrick scaling requirement. (See footnote 1.) After the introduction of the 
Chern-Simons (CS) term (and of Aq) however, such solutions were not found. This is because the appearance 
of the CS term in the action results in the effective presence of a higher order term in the static energy density 
functional, as explained in Section 1. This effective higher order CS term, which scales as then destroys 
the previous Derrick scaling balance. This is why we have proposed models in d = 2n + 1 for n > 3 only, 
since in d = 5 the only choice for the ( "magnetic" ) YM sector is the p = 1 YM term. How this obstacle is 
circumvented in the case of LT dyonic instanton |14j (and all higher analogues in d = Ap + 1 dimensional 
spacetime) can be seen in Appendix C. 

As a final comment, we note that the solutions discussed here are exclusively spherically symmetric. One 
would expect that there exist dyons of these systems subject to less stringent symmetries, though technically 
the numerical construction of these putative dyons would be a formidable task. 
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monopoles on MP . But the existence of a scalar valued global electric flux relies on the definition of a ' Hooft electromagnetic, 
something that is problematic in flat spacetime dimensions higher than 4. 

5 Though in that case one can employ instead any one of the excited Georgi-Glashow models, or any superposition thereof, 
which descend from higher than 4 dimensions p]. 
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A Static limit of the C hern- Simons density in 2n + 1 dimensions 

To illustrate our statement, namely the result p,15[) . we carry out the relevant computations explicitly in 

)(2n- 
'cs 



3, 5 and 7 dimensional spacetimes. The CS densities in question, for n = 1, 2 and 3, are 



(3) _ 

cs — 



£Au,Tr A 



^cs — ^A/ii/po-Tr -\ A 



F _ Z. A A 

1 fiv 2 A* ^ 



(A.l) 



(7) 



F F F p F A A F A F A 

± pv- pa- th ^- pv- pa- rL T- rL K ± pv--p- UT^K 



"^F^ v ApA a A T A K — -^A^A v ApA a A T A K 



which are by construction gauge variant. Their variations however are by contrast gauge covariant and are 
easily expressed in general 2n + 1 case as 



- q(2tI+1) _ / , -I \ p A|il|*2A«3M4-M2n-l/*2B F F F 

°'4a 1 'CS ~~ r piP2 T P3Pi ■ ■ ■ F P2n~\P2~ 



(A.2) 



The statement here is the following: for static fields, the CS density on d = 2n + 1 dimensional 

spacetime reduces to iS.15\) , namely to 



n 



(5) 
CS 



fig s n+1) = V ■ Q. {2n+1) + (n + l)Tr F A F A • • • A F X, n times . 
To illustrate this, we give the explicit total divergence expressions for n — 1, 2 and 3, (|A.3[) are, 

fig] = 2 £ij Tr | A F, 3 j + V • O (3) 

= 3 £ ijH Tr j A) F y F fei J + V • n< 5) 

= 4 Eijklmn | ^0 -Fij i^fei ^mn | + V • , 

with ^ 2 ™ +1 > = Of n+1) given by 

-2 e^TrjAo Ajj, 

-2 e ijfcI Tr j A [(A, F u + F u A,) - Aj A k A{\ 



(A.3) 



(A.4) 



"cs 



n 



(3) _ 







(5) 



(A.5) 







(7) 



-2e i7 - feJ TWA 



(-Fy -Fw + Fij A n Ff~i + A n Fij F k t) 



- - (A k A t A n Fij + F^ A k A t A n ) 

O 

2 4 
--(A k At F^ A n + A n F^ A k At) + -Ai Aj A k At A % 



D 



B Spherically symmetric Ansatz for SO(D + 2) YM on R 

The static spherical symmetric Ansatz for the SO(6) Yang-Mills (YM) system on R 6 , used in [T5], is extended 
to the SO(D + 2) YM on H D , i.e. in D + 1 dimensional spacetime with D — In. Since D is even, we can 
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take the YM connection to take its values in one or other chiral representation of SO±(D + 2), such that 
our spherically symmetric Ansatz is expressed in terms of the representation matrices 



-(±) 

J a/3 



i ( i±r 2n + . 



[T a> T ] 



a,p = l,2,...,2n+2 , 



T a = (Ti, Tm), with the index M = (2n + 1, 2n + 2), being the gamma matrices in 2n - 
T2n+3, the corresponding chiral matrix. 

Our spherically symmetric Ansatz for the YM connection = (Ao,Ai) is 



(B.l) 

2 dimensions and 



A Q = 

Ai = 



-( £ x) M ^-4m-X 2 " + - 1 



J ij 3 



2n+l,2n+2 

< / )M \ M 



s (±) + 



-A 



2n+3 ~. £<(±) 



« ^ J 2n+l,2n+2 



(B.2) 



(B.3) 



in which the summed over indices M, N = 2n + 1, 2n + 2 run over two values such that we can label the 
functions (</> M , </> 2n+3 ) = 0, (x M ,X 2 ™ +3 ) = X and (Af ,A?"+ 3 ) = A r like three isotriplets $, x and l r , all 
depending on the 2n dimensional spacelike radial variable r. e is the two dimensional Levi-Civita symbol. 

The parametrisation used in the Ansatz (|B.2j) - (|B.3j) results in a gauge covariant expression for the YM 
curvature F^ v = (Fij,F iQ ) 



1 



D r 



Ft 



ai 



b ^+3 + 1 ^|2 _ j 
x 2n+3 



tMi v(±) 
5 ^[^jJM 



2n+3 « 



(±) 



in which we have used the notation 

D r <j) a = d r <j) a + e abc A\ (/> c 



D rX a = d rX a + e abc A b rX 



(B.4) 



(B.5) 



(B.6) 



as the 5*0(3) covariant derivatives of the two triplets = (j) a — ((f> M , (j) 2n+3 ), x = X a — (x M :X 2 " +3 )j with 
respect to the 5*0(3) gauge connection A r = A". This gauge connection (in one dimension, r) is of course 
trivial and hence does not appear in the field equations, but it serves to give the constraint equation. 



C Dyonic instantons in d = 4p + 1 dimensional spacetimes 

The dyonic instanton is a (static) solution to the Yang-Mills-Higgs model in 4 + 1 dimensions. It shares with 
the Julia-Zee (JZ) dyon, the presence of a Higgs term in the action, which is to be identified with the electric 
non-Abelian connection A leading to the definition of the electrix flux. This is in direct contrast with the 
models proposed here, which feature no Higgs field and Aq appears through the Chcrn-Simons terms in the 
action. 

In contrast to the JZ dyon however, static energy density functional of the Lambert-Tong (LT) dyonic 
instanton does not depend on the Higgs field and consists solely of the usual Yang-Mills action density. 
More specifically, the YM field is the (anti-)self-dual instanton on R 4 . As such the Higgs field plays an 
auxiliary role in the description of the LT dyonic instanton, and is needed only for the definition of the 
global electric flux. This makes it possible to generalise the construction of dyonic instantons to 4p + 1 
dimensional spacetimes, albeit in a rather limited context. This is our aim in the present Appendix. 

We limit our considerations here to the recovery of the 4 + 1 dimensional LT dyonic instanton, and 
proceed to its 8 + 1 dimensional version. The Ap + 1 dimensional case then follows systematically. The 
Lagrangian Ld+i for D = 4 is 

£4+i = ^ yy{i;,l" } + -T^A*^*} , (C.l) 
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yielding the following gauge field equations 

ZV^o-[$,A)$] - (C.2) 
DpF^-DoFov-faD,,*] = 0, (C.3) 

and the Higgs field equation 

D 2 § - Dl<f> = (C.4) 

having used the Minkowskian metric rjij = diag(+, — , — , — , — ), with 7700 = 1 and rj^ — —5^, i — 0, /i. 
Eqn. (|C.2[) is the Gauss Law equation, and (IC3|) the Ampere equations. The Hamiltonian is 



iTrH 
4 \2 

i_ r 1 



^4+1 = 7 Tr\-F^ + Fl + D^ 2 + D ^ 2 \ (C.5) 



4 Tr + (F 0M - + D & + 2D^F A . (C.6) 

Following [TJ], we have expressed the second line of (|C.6p such, that after substituting the Gauss Law 
equation (|C.2[) and discarding the surface term arising from 

a M Tr{$F 0A1 } , 

(|C6|) takes the form 



H i+1 = \^i\Fl v + {F 0li -D^f + D ^\ . (C.7S 



We now restrict our attention to static solutions, and, choose 

A Q = $ = , F 0Al - £> M $ , (C.8) 

as a result of which the static limit of the Hamiltonian (|C.7j) just reduces to the action density of the 1-YM 
system, i.e. only the first term in (|C.7j) . whose action is minimised absolutely by the 1-BPST solution [TO] . 
Using (|C.8I) . one can readily verify that the Ampere Law equation (|C.3I) reduces to 

= 

which is solved by the selfdual BPST fields, and the Gauss Law equation (|C.2I) and Higgs equation (|C.4I) 
become identical, reducing to 

2V?„$ = 0. (C.9) 



There remains to find a regular solution to (|C.9[) . which we delay till after treating the 8+1 dimensional 
case, next. 

The Lagrangian £d+i in the D = 8 case is 

£s + i = -^Tr{F ijkl F^ kl } + JL T r{F ufc F^ fc } (CIO) 
in which we have used the notation 

Fijki — {Fi[j,F kl }} = {Fij,F k i} + {F ik ,Fij} + {Fii : F jk } (Gil) 
F ijk = {%, D k] <£} = {F&D**} + {F jk ,Di<f>} + {F ki ,Dj$} . (C.12) 

The choice of the Higgs kinetic term Ti:{Fuj, D k ^} 2 here, instead of the usual one TrZ?i$ 2 , is made so that 
the static Euler-Lagrange equation for A becomes identical to the corresponding Higgs equation after the 
identification A — $ is made. 
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Using the same Minkowskian metric as above (but in 8 + 1 dimensions now) we have the following 
Yang-Mills equations 

D li {F p(T ,F 0fipa } + 2D fl {D p ^,F ^} + [^,{F pa ,F 0pa }} =0 (C.13) 

Dfj,{Fp a , F avpa } + 2£» M {£» P $, F^p] - 2D„{F 0p , F ^ p } - £> M {D $, F „ v } 

-D {F pa , F Qupr7 } - 2D {D p $, F^p] 

-[$,{F p<7> F vp<T }}-2[$,{F 0p ,F 0vp }} = (C.14) 

and the Higgs equation 

D P {F pa ,F ppa } + 2D p {F 0v , F 0pv } - D {F pv , F 0pv } = . (C.15) 

Eqn. (IC.13I) is the Gauss Law equation and (|C.14p the Ampere equations. 
The Hamiltonian of (jC.101) is 

^ 8 +! = 77 Tr 1 1 F tv P a + F 0uvp + F'l v + 2Fq 1 (C.16) 



|^ ^ pupa ' ± Opvp ' ± ^ivp ' Qpi> 

-A- 

4! 1 4 



Tr l ~ F 1vp° + (-^om^p - ^p) 2 + 2 F 0p,u + 2F Qpvp F pvp \ (C.17) 



(C.18) 



where again, as in the D = 4 case above, we have expressed the second line of (IC.17|) such, that after 
substituting the Gauss Law equation (jC.13l) and discarding the surface term arising from 

d p Tr<P{F^,F ppu } 

(|UT7|) takes the form 

H 8+ i-^Tr ( 

Again, we restrict our attention to static solutions, and as in the p = 1 case above we choose 

A = <S>=^ D <S> = , F 0/1 = F tlvp =F 0piVp , F Om „=0, (C.19) 

as a result of which the static limit of the Hamiltonian (jC.18l) just reduces to the action density of the 2-YM 
system, i.e. only the first term in (|C.18|) . whose action is minimised absolutely by the 2-BPST solution 
Using (IC.19[) . one can readily verify that the Ampere Law equation (|C.14[) reduces to 

F) p {F p(71 F pupa } 

which is solved by the selfdual p = 2 BPST fields, and the Gauss Law equation (jC.13l) and Higgs equation 
(|C.15P become identical, reducing to 

D^{F p(7 ,{F [p(T ,D^}} = 0. (C.20) 

Up to this point, the D = 4 and the D — 8 cases are on the same footing, and hence also is the generic 
D = Ap case. The remaining task in all these cases is the construction of a regular solution to the Higgs 
equations (|C.9|) , (|C.20|) , etc. , such that the magnitude of the Higgs field tends to a constant at infinity. This 
is to guarantee the convergence of the surface integral 



/ TY{<S>F p0 }dS», (C.21) 

Js - 1 



for the electric flux. The integral (|C21|) will be convergent provided that the Higgs field $ = A decays 
fast enough at infinity, and, is regular at the origin. Such a solution was found p3] for D = 4 (p = 1) but 
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unfortunately, we have not succeeded to do this for the new, p > 2 cases. Let us discuss this question a little 
further. 

In the D = 4 (p = 1) case, the Higgs equation (|C.9[) is solved most conveniently in the background of the 
BPST field in the 't Hooft singular gauge 

4 ( f)= l^(r) E £v u;(r) = _^!_l! j (c.22) 

by positing the following Ansatz for the Higgs field 

$ = h(r) £ 34 , (C23) 
and noticing |^| that the inverse of the function multiplying £ M „ in (|C.22I) , namely that 

Hr) = r (C.24) 
1 — w(r) 

solves (IC.9|) . (In (|C.22|) - (|C.23|) . £ M „ are the chiral representations of the algebra of 5(0(4).) 

The situation with solving (|C20[) in the background of (|C.22[) (now interpreting Y,^ there as the chiral 
representations of the algebra of SO(8)) is a very difficult task, which we have not succeeded in. Even 
finding an Ansatz like (|C.24|) for $, which yields a single ordinary differenial equation (ODE) for h(r) is 
problematic. 

Recognising that the Higgs field (or Aq) does not necessarily have to take its values inside the algebra of 
the "magnetic" field, we found an Ansatz for the Higgs field that leads to a single ODE. This is 

$ = h(r) S 56 78 , S 56 78 = {S 5 [6, S 78] } ee {S 56 , S 78 } + cycl. (678) , (C.25) 

for which the equation (|C.20[) reduces to one single ordinary differential equation for the function h(r). This 
is 

- 16A 2 (3A 2 + 2r 2 )h + 3r(A 2 + r 2 ) ((7A 2 - r 2 )ti + r(X 2 + r 2 )h") = 0. (C.26) 
The asymptotic solutions of (|C26|) are found. As r — > 0, h(r) behaves as 

h(r) = g 2 r 2 gfr 4 + |§_r 6 + 0(r 8 ) (C.27) 

(with a free parameter (72), while it decays at infinity as 

W=ft + ^ + ^ + 0(l/r») (C28) 

go being another free parameter. 

Unfortunately, there exist no solutions h of (|C.26|) that satisfy the two asymptotic values in (IC.27I) and 
(|C.28|) . simultaneousljQ. This can be shown as follows. Let gi > 0; then according to (IC.27|) h is increasing 
for small r. Choosing go > in (jC.28l) for large r, h' < in this region, so that 52 > go- This means 
that h must have a positive valued maximum, h maK > in between the two asymptotic regions. At this 
maximum h' max = and h'^ax < 0. Substituting /i max > 0, h' max — and h'^ lax < in (|C.26[) . one sees that 
the nonvanishing terms there are both negative, and hence (|C.26|) cannot be satified. (This argument can 
be reversed starting with g 2 < 0.) 

So even with the special Ansatz (|C.25[) . there exist no solutions supporting a nontrivial Higgs field for 
the p > 2 case presented above. Since the Higgs field is necessary for the definition of a (finite) global electric 
charge, this means that we cannot find a natural generalization of the LT dyonic instantons in d — 4p+ 1 > 5 
dimensional Minkowski space. 



"Likewise, identifying the prepotential function in the Jackiw-Nohl-Rebbi 22 Ansatz with the inverse of the corresponding 
Higgs field prepotential, yields the multi-dyonic instantons 1231 . 
7 We thank J. Burzlaff for help on clarifying this point. 
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